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Abstract. Due to the presence of light pions in the theory, lattice QCD at finite densities 
suffers from issues with noise in both grand canonical and canonical formulations. We study 
two different formulations of the Nambu-Jona-Lasinio model reduced to 2+1 dimensions at 
large N, where iV is the number of flavors. At finite chemical potential one formulation has 
a severe sign problem and a fermion correlator which displays a broad probability distribution 
with small mean. In the other we find no sign problem and a distribution amenable to the 
cumulant expansion techniques developed in Ref . [TJ [21 E] - 



1. Introduction 

While lattice QCD is currently the only tool for studying QCD from first principles in the 
non-perturbative regime, the study of QCD at finite densities on the lattice has been greatly 
hindered by two related problems. In the first, lattice QCD at fixed baryon chemical potential 
suffers from what is known as a sign problem, which arises due to the complexity of the fermion 
determinant one wishes to use as a probability measure for Monte Carlo calculations. A possible 
solution is to employ phase reweighting, in which the phase of the determinant is absorbed into 
the observable and the magnitude of the determinant is used as a probability measure; however, 
the expectation value of the phase using the same ensemble must also be calculated in order 
to reconstruct the original observable. One finds that for /i above some critical value, this 
expectation value vanishes exponentially and becomes swamped by statistical noise. The second 
problem is encountered in a canonical formulation, in which a fixed number of quark sources 
and sinks are separated by a long Euclidean time to project out the ground state. Again, the 
signal-to-noise ratio is found to vanish exponentially with time. 

A notable property of these two problems is that they are both related to the same dynamical 
property of QCD; namely, they both arise due to the presence of light pions and thus are 
intimately related to the phenomenon of chiral symmetry breaking. In the case of finite chemical 
potential, this can be seen by noting that for two degenerate flavors of quark the magnitude of the 
fermion determinant corresponds to chemical potential for isospin. Thus, at low temperatures 
for fi > m-n/2 pion condensation would tend to occur, and the role of the phase must be to 
cancel this pion condensation [31 [5] . 

For the other case where finite baryon density is achieved using correlation functions of quark 
propagators, the role of the pion in the signal-to-noise problem may be noted by considering the 



variance of an operator used to create nucleons. While at large Euclidean times the correlator for 
three appropriately contracted quark propagators projects out the ground state of the nucleon, 
Mat, the variance will be composed of three quark and three anti-quark propagators, projecting 
out the pion mass at large times [6J [7]. The signal-to- noise ratio will thus be exponentially 
vanishing with the factor (M^y — |m, 7r )r, where r is the Euclidean time, for late Euclidean 
times. 

Both of the noise problems presented above can be shown to arise from probability 
distributions which are broad, symmetric, and centered about zero. For the nucleon correlator, 
one may see this by looking at higher moments of the distribution function. Even moments 
consist of equal numbers of quark and antiquark propagators which may pair up to form pions, 
while odd moments have three unpaired quark propagators which will form a baryon [8] . Thus, 
because odd moments will always contain a baryon mass in the exponential, they will be 
suppressed relative to the even moments, giving a symmetric distribution with exponentially 
small mean. A similar argument may be made for the expectation value of the phase at 
finite chemical potential, where all odd moments are equivalent to the mean, and therefore 
exponentially small for fi > m 7r /2, while all even moments give exactly 1 (for a discussion of the 
distribution of the phase, see [9]). 

In [H El [3] it was shown that another type of noise problem exists in which the distribution 
displays a long tail, resulting in a sample mean which may be far from the true mean with 
deceptively small error bars. This problem is often referred to as an overlap problem, and 
has been known to occur for observables when reweighting is performed. This was discussed 
in a canonical formulation using lattice calculations of unitary fermions, a strongly-interacting 
nonrelativistic system in which the fermion determinant is real and positive for even numbers of 
fermion flavors. The probability distributions for fermion correlators were shown to be nearly 
log-normal, and a possible solution to the overlap problem was postulated in the form of a 
truncated cumulant expansion for the correlator, 
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where 0(4>) is some operator of interest on a background field (ft, n n (lnO) is the nth cumulant 
of the logarithm of the operator, and N K is some suitable truncation order, determined by 
examining the convergence of the series. For a distribution that is sufficiently close to log- 
normal, this series may converge at a low enough order to be useful in practical applications. 
Recently it has also been reported that nearly log-normal distributions are also seen in lattice 
QCD calculations of correlation functions for intermediate times [10]. 
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2. Two formulations for a large- N NJL model in 2 + 1-dimensions 

To understand the origin of these two types of noise problem in QCD and to further explore the 
link between signal-to-noise and the presence of light pions, we turn now to a QCD-like model 
which displays chiral symmetry breaking without the added complication of confinement. For 
more details behind these calculations, see [11]. The model we will consider is the NJL model 
in 2 + 1-dimensions for large number of fermion flavors, N, 

jO = NU a (0-m)i/j a -^[^ a i/j a ) 2 + ^ a i^ a ) 2 ]) , (2) 

where a, b... are a flavor index summed over 1, ...,N, indices, are summed over 3D 

coordinates 1,2,3, while Greek indices fi, v... are summed over 4D coordinates 1, ...,4. The 
gamma matrices are the usual 4x4 matrices used in 4D, and so the Lagrangian represents 2N 



flavors of 3D Dirac fermions. In the limit m — > this theory has a U(l) chiral symmetry in 4D, 
which becomes a flavor symmetry in 3D; this "chiral" symmetry is spontaneously broken as in 
3 + 1-dimensions, giving rise to a Goldstone boson. 

To study this model numerically one may introduce auxiliary fields to generate a four- 
fermion interaction, which can then be integrated over using Monte Carlo methods. We find two 
physically equivalent methods for doing so. The first, which we will call the a/n formulation, is 
the conventional method which introduces scalar fields that are singlet under the SU(2N) flavor 
symmetry, 



C = N f (a 2 + 7T 2 ) + [0 - m + a + ivr 75 ] ^ 
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The second formulation, which we will call the A/V formulation, follows if one performs a Fierz 
rearrangement of the four-fermion interaction in eq. ^ before introducing auxiliary fields. This 
results in the introduction of N x N matrix valued vector and axial vector auxiliary fields V 
and A, and the equivalent theory 

£ = N Tr (VM + A M + ^a[0-m + if + 4 75 ] ab fy\ . (4) 

To determine whether these two formulations will cause a Monte Carlo sign problem at non- 
zero chemical potential, we should investigate the positivity of the fermion determinants. For the 
ct/tt formulation, it is possible to define a real symmetric charge conjugation matrix C satisfying 
C 2 = 1, C'y^C = 7* for i = 1,2,3 and C75C = —75 . Then the fermion operator for a single 
flavor in the grand canonical formulation satisfies D* = CDC, where D = 

and complex eigenvalues of D must come in conjugate pairs. Thus (det D) N is real, and positive 
for even N. This implies that there is no sign problem at finite density |12| . 

In the A/V formulation, the fermion matrix at finite chemical potential is given by D(fi) = 
(0 — m + iY + ^(75 + /U71) , which is similar in structure to the QCD Dirac matrix with nonzero 
/i, and its determinant is similarly complex. In fact, as in QCD, the magnitude of the fermion 
determinant for two degenerate families, | det D(n)\ 2 , corresponds to isospin chemical potential, 
so that for \i > m n /2, the phase is responsible for eliminating pion condensation in the ground 
state. 

Chiral symmetry breaking is simple to see in the a /it formulation, where the large N 
expansion is equivalent to the semiclassical expansion. Upon integrating out the fermion fields 
we have the following effective action: 
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Performing a mean field calculation leads to (using dimensional regularization and the MS 
subtraction scheme), 

(a) = {-\l + y/l + Am/ f + 2m/ f] , M = {a) - m = f - \l + y/l + Am/f] , (6) 
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where we have defined (a) = f to be the chiral symmetry breaking minimum when m = 0, 
and M is the constituent fermion mass. We may also find the a and tt dispersion relations 
D a (k), D n (k) by expanding the effective action to second order about our mean field solution. 

In the A/V formulation we cannot use mean field theory; instead, to leading order in 
1/N we may find the fermion propagator by solving the Schwinger-Dyson equation and 
finding the nonzero fermion mass, Eq. [6} Furthermore, one may derive the a and it 
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Figure 1. Graphical representation of the integral 
equation for the four-point correlator in the A/V 
formulation for an incoming fermion-antifermion 
pair of one flavor and an outgoing pair of another. 
Dirac indices are labeled. 

meson propagators by solving the matrix equation shown in Fig. [TJ One finds Aiij-ki(k) = 
— (Sij5ki/D (T (k) + (i , f5)ij(ij5)kl/D n (k)). Thus we see that an interaction via i-channel exchange 
of and mesons is exactly equivalent to a single meson in the s or u-channel in the a /n 
formulation, corresponding to a valence fermion/antifermion pair or two valence fermions or 
antifermions, respectively. 

3. Probability distributions for fermion correlators 

We now proceed to calculate the probability distribution for a fermion propagator in both 
formulations. If X[<f>] is a functional of a stochastic field <f> corresponding to an observable, 
we define the normalized probability density function for X to be the path integral V(x) = 
Af J[d<f>] e~ s W 6(X[<I>] - x), where we assume S[</>] is real. This quantity is difficult to analyze 
field theoretically due to the singular nature of the delta function, so we perform a Fourier 
transform and work with the characteristic function, 
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where we have generalized the cumulant generating function W(s,s) = — ln<3?(s,,s) for a 
complex observable through a double expansion in both s and s. From this formulation we 
see that we may compute cumulants using connected Feynman diagrams of the modified action 
S[</>]-i(8X[4] + 8X[<f>]). 

For our observables we will use 

a /it formulation, 

(8a) 

A/V formulation, 

(8b) 

where V is some Dirac matrix of our choosing, and we consider the logarithm of the propagator 
in the a/n case to simplify the calculation later on. Measuring the expectation value of this 
correlator is a procedure for determining the mass mj of the lightest fermion state allowed for 
a given T using limT->oo — y ln(Xr) = limT->oc — y ln(e yr ) = m/. 

The calculation of the variance for the A/V case involves attaching fermion propagators at 
zero spatial momentum to the legs in the first diagram in Fig. [TJ For T = 1 and near the chiral 
limit the pion is the lightest state which can propagate through such a graph, so that for late 
Euclidean time we find, K\ \ c± — 2 J[ /M e~ m " wT . Higher cumulants can be computed using the 

equivalent a /it diagrams. Generically, a graph contributing to K m ,n scales as JV m+n_1 , since 
we need a minimum of (m + n — 1) mesons to make a connected graph. Furthermore, it is 
straightforward to verify that the minimum mass state that can possibly propagate in a graph 
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for K mn with m > n consists of (to — n) fermions with mass M and n pions. Therefore we 
expect these cumulant to scale as 
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(m > n) . (9) 



The above scaling implies that the distribution for the real part of the fermion propagator near 
the chiral limit becomes highly symmetric about zero at late time. That is because odd moments 
(for which (to — n) ^ 0) are seen to fall off much more quickly than even moments. This is 
consistent with the Lepage-Savage picture for baryon propagator distributions in QCD. 

For the a/ir formulation the cumulants are given by the connected graphs derived from the 
action Sy = NS[a,ir] — isYr[a,ir], where S[a, tt] is given in eq. §5§. At leading order in a 1/N 
expansion, K n will be given by the sum of all tree level diagrams composed of n external lines 
and any number of vertices which arise from internal fermion loops. Given our definition of Yp, 
we find that all external propagators are spatially homogeneous, and, for very large Euclidean 
time, temporally homogeneous up to edge effects near the source and sink. Because we only 
consider tree graphs, all internal meson propagators will similarly be spatially and temporally 
homogeneous; thus, we may once again use mean field theory. Solving for the mean field leads 
to the following cumulants for the logarithm of the fermion correlator: 

„ = In,- MT, K „ i2 = -<|^|i!(2M-/)TC»-', { * ^ _ /)2 (1Q) 

where z gives the overlap of our operator with the ground state, and we have approximated 
the integral over Euclidean time by a step function for —T/2 < r < T/2. Note in the limit 
N,T,V — > oo, yjj = finite, the cumulants K n vanish for n > 3 and Yp assumes a normal 
distribution, giving a log-normal distribution for the correlator. With the variance for the log 
of the correlator growing linearly with time, the distribution for the correlator will eventually 
become heavy-tailed, with a skewness growing exponentially with a, though such long times are 
unlikely to be reached in practical numerical simulations. 



4. Discussion 

We found that in the QCD-like U A/V" formulation of the NJL model, the fermion determinant 
was complex and a Splittorff-Verbaarschot argument [U [5] could be made to show that the 
phase of the fermion determinant had to fluctuate wildly for /j, > m n /2. When looking at 
fermion correlators, the distribution evolved to be symmetric about an exponentially small mean 
relative to its width, implying a severe signal/noise ratio when sampling the correlator using 
Monte Carlo methods. In contrast the "cr/V" formulation with even N has no sign problem at 
nonzero fi, and the correlator distribution was, in a certain limit, log-normal and heavy-tailed. 
This overlap problem would pose challenges to Monte Carlo sampling if the tail became too 
long, but this sort of problem seems to be less severe than the exponential fall-off of signal/noise 
seen in the A/V formulation as seen with the cumulant expansion analysis of Ref. pi 1131 [14]. 

Our analysis should make it clear that the sign problem encountered in QCD at nonzero 
chemical potential is not a fermion problem, but instead a consequence of interactions. In 
particular, if the particles being studied can exist in a tightly bound state of valence fermions, 
there is going to be a sign problem. We believe that inventing a way to introduce the pion into 
QCD as a fundamental field could be an important step toward solving the QCD sign problem 
and beginning to study the properties of ordinary and dense matter from first principles. 
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